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Abstract 

A new class of exact solutions of Einstein's equations is derived pro- 
viding a physically plausible hydrodynamic description of cosniological 
matter in the radiation era of the Universe expansion. The solutions are 
characterized by the LTB metric with dissipative fluid sources, subjected 
to the following conditions: (i) the nonequilibrium state variables satisfy 
the equations of state of a mixture of rclativistic and non-relativistic ideal 
gases, where the internal energy of the latter has been neglected, (ii) the 
particle numbers of the mixture components are independently conserved, 
(iii) the viscous stress is consistent with the transport equation and entropy 
balance law of extended irreversible thermodynamics with the coefficient of 
shear viscosity provided by kinetic theory. The Jeans mass at decoupling 
is of the same order of magnitude as that of baryon dominated models (i.e. 
Mj w IQISa/q). 



1 Introduction 

The standard approaches to the radiative era of cosmic evolution (the period 
from the end of cosmic nucleosynthesis to the decoupling of matter and ra- 
diation) resort to a FLRW space-time background while the sources of the 
gravitational field are described either by equilibrium kinetic theory Q, gauge 
invariant perturbations or some hydrodynamical model |3[, which, in gen- 
eral, fail to incorporate a physically plausible description of the matter-radiation 
interaction since they assume thermodynamical equilibrium throughout. Here 
we consider the aforesaid sources in the temperature range lO^K < T < lO^K, 
as a nonequlibrium mixture of a non-relativistic fluid (matter) and a extreme 
relativistic one (radiation); their mutual interaction is modeled by a dissipa- 
tive shear-stress tensor, and the background space-time is described by the 
Lemaitre-Tolman-Bondi (LTB for short) metric @. (The interest for inhomo- 
geneous metrics has revived in the wake of the work of Mustapha et. al. who 
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show that given any spherically symmetric geometry and any set of observa- 
tions, evolution functions for the observed sources can be found that will make 
the model in general compatible with observation). The stress-energy tensor 
reads 



nab 



where is assumed that the matter fluid is just "dust" and therefore the hydro- 
static pressure is furnished by the radiation only, i.e. 

In its turn the shear dissipative pressure tensor is governed by 



+ 2ri Gab = 0, with r/, 



' (r) 



this expression is compatible with relativistic causality and stability, and is sup- 
ported by kinetic theory, statistical fluctuation theory, and information theory 
H. The entropy per particle takes the form 



nl 



(snn'^);a > 0, 



with 



a 



and r is the relaxation time for the shear-stress. Further, since the number 
of particles of each fluid is independently conserved, we have (n^'^^n");^ = 



2 Exact Solution and Density Contrasts 

The LTB metric can be written as 



ds^ = -c^dt^ + - -dr'^+Y'^ ide'^ + sm^{e)d(t)'^] , Y = Y{t,r), F = F{r). 

\ — F ^ ' 

From this metric element the expansion and shear read 

o"=-|— -— I, 0"";, = diag [0, -2a, a, a] , 



y 2Y 
= 77 + 1^ 



3 \ y Y' 



respectively; and the most general form for the shear-stress reduces to H'^ ^ = 
diag [0, —2P, P, P], with P = P{t, r). The non-trivial Einstein's field equations 
are 

Y + i7c2^ 



Kp 



Y2Y' 



YiY'^ + Fc^] +2y2y 



3y2y/ 



U2G% + G\) 
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kP 



Y 

W' 



Y [Y^ + Fc^] + 2Y^Y 



Y3 



where the prime and overdot mean partial derivatives with respect to r and t, 
respectively. In the case of flat spatial sections {F = 0), integration of these 
equations yields the following exact solution |0] 



-^{t-ti) = ,/V+l{y-2e)-^T+~e{l-2e), 



where we have used the definitions 



Y? ' 



W_ 

M' 



y = 



Y 
Y, 



with 
M = 



Pi Yi Yidr, 



(m) 2 

pl = mc 



W 



Jp^^Y^Y^dr, pf-^ ^Snt'k.T,, 



so that p["^^ , p[^^ define the initial densities of the non-relativistic and relativis- 
tic components of the mixture, respectively. 

At this point it is expedient to introduce the density contrast parameters 



(m) 



1 + A 



(m) 



Pi = \Pi 

as well as the ancillary functions 

y/y T (1-r) 



1 + A 



(r) 



Y,'/Y, 3(1 + A|'-')' 

the latter are linked to the density contrasts by 



$ = 1 + 



(1 - 4r) 



3(1 + A 



r = 1 + 3^A^'"' + 35 Af' , 



AAA'-^'' + (1 - AB) Af ' 



1 + a; 

where the quantities A and B are known functions of y. Given a set of initial 
conditions specified by e, A^™', A^"^', on some initial hypersurface, the forms 
of the state and geometric variables as functions of y and the chosen initial 
conditions are fully determined. However, for the solutions to be physically 
meaningful they must comply with the following set of physical restrictions, 

. 15^5 f<^\ 



16r V^J 



f > 0, 



2aT [P. 



s 3^<^ „(0 

ri 



1 + 



3p 



(r) 



<0, 



r < 3/G before decoupling, r > 3/0 after decoupling. 
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3 Initial Perturbations and Jeans Mass 



Assuming that the perturbations on the intial hypersurface are non-adiabatic 



(i.e. ^ = |A['^ - Af"'' 7^ 0) and that Ti « 10^ Kelvin, the set of values asso- 
ciated to the matter-radiation decoupling is obtained from solving numerically 
the equation [0] 



(m) 



r,-4xl03 = i^^(y,,AW,Ar: 



(T 



which results in y/j ~ 10 



,2.4 



The temperature anisotropy of the CMR ([^T/T]^ ~ 10 ^) constrains the 



maximal values of A,-"', A-*^' to about 10 ^ and the corresponding variation 



range of |A 



is lA.^'^l < 10" 



and so compatibility with acceptable values of 



and the CMR anisotropy implies |A, 



(«)i 



iArp<iAn 



For non-adiabatic perturbations the Jeans mass is given by 



47r 



-mn 



(m) 



k{p + p) 



1 3/2 



An c^xor 



1/2 



3/2 



where 



1 + 



3p(™) 



Assuming y = yn ^ 10^'^, e l/xi ~ 10^ and pf^ ^ o.s'^t ^ 7.5x10^ ergs/cm^, 
yields Mj ~ lO^^gm, or approximately 10^^ solar masses. This value coincides 
with the Jeans mass obtained for baryon dominated perturbative models as 
decoupling is approached in the radiative era. 



4 Conclusions and Outlook 

The model presented here is a self consistent hydro dynamical approach to 
matter-radiation mixtures that: (a) is based on inhomogeneous exact solutions 
of Einstein's field equations, and (b) is thermodynamically consistent. It has, 
however, the limitations of not incorporating heat currents (the LTB models 
do not allow for thisQ) nor bulk dissipative stresses. Nonetheless, we believe 
this model may become a useful theoretical tool in the study of cosmological 
matter sources, providing a needed alternative and complement to the usual 
perturbative or numerical approaches. 

The solutions have an enormous potential as models in applications of as- 
trophysical and cosmological interest: 
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Structure formation in the acoustic phase. There is a large body of hterature 
on the study of acoustic perturbations in relation to the Jeans mass of sur- 
viving cosmological condensations. Equations of state of the type "dust plus 
radiation" are oftenly suggested in this context Q. Since practically all work 
on this topic has been carried on with perturbations on a FLRW background, 
the exact solutions presented here may be viewed as an alternative treatment 
for this problem. 

Comparison with perturbation theory. Our model is based on exact solutions 
of Einstein's field equations, but their initial conditions and evolution can be 
adapted for a description of "exact spherical perturbations" on a FLRW back- 
ground. It would be extremely interesting, not only to compare the results of 
this approach with those of a perturbative treatment, but to provide a physically 
plausible theoretical framework to examine carefully how much information is 
lost in the non-linear regime that falls beyond the scope of linear perturba- 
tions. We have considered only the case -F = 0, thus restricting the evolution 
to the "growing mode" since all fluid layers expand monotonously. The study 
of the more general case, where F{r) is an arbitrary function that could change 
sign, would allow a comparison with perturbations that include also a "decay- 
ing mode" related to condensation and collapse of cosmological inhomogeneities. 

Inhomogeneity and irreversibility in primordial density perturbations. The ini- 

tial conditions of the models with 7^ are set for a hypersurface with 
temperature Tj ~ lO^K. These initial conditions can be considered the end 
product of processes characteristic of previous cosmological history, and so the 
estimated value |A|^^| ~ 10^^, related to the spatial variation of photon entropy 
fluctuations, can be used as a constraint on the effects of inhomogeneity on pri- 
mordial entropy fluctuations that might be predicted by inflationary models at 
earlier cosmic time. Also, the deviation from equilibrium in the initial hyper- 
surface (proportional to [A-'']^ « l^l**^!) might be helpful to understand the 
irreversibility associated with the physical processes involved in the generation 
of primordial perturbations . 
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